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Abstract 

We extend our previous work on the cosmology of Coleman-de Luccia bubble colli- 
sions. Within a set of approximations we calculate the effects on the cosmic microwave 
background (CMB) as seen from inside a bubble which has undergone such a collision. 
We hnd that the effects are always qualitatively similar — an anisotropy that depends 
only on the angle to the collision direction — but can produce a cold or hot spot of 
varying size, as well as power asymmetries along the axis determined by the collision. 
With other parameters held fixed the effects weaken as the amount of inflation which 
took place inside our bubble grows, but generically survive order 10 efolds past what 
is required to solve the horizon and flatness problems. In some regions of parameter 
space the effects can survive arbitrarily long inflation. 
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"Through telescopes men of science constantly search the infinitesimal corners of our 
solar system seeking new discoveries, hoping to better understand the laws of the Universe. " 

- When Worlds Collide (The Movie) 

1 Introduction 

Perhaps the greatest challenge for string theory is to find a testable prediction which can 
differentiate it from other potential theories of quantum gravity. This has proven to be 
extremely difficult, because the scale at which string theory effects become important is be- 
lieved to be far out of reach of the energies accessible to any earth-based particle accelerator. 
However cosmology provides another opportunity: the early universe was arbitrarily hot and 
dense, and remnants of the physics of that time remain in the large-scale structure of the 
universe today. 

It has become increasingly clear in recent years that one of the characteristic features 
of string theory is its so-called landscape [H El [3]. The landscape is a potential energy 
function of the would-be moduli scalar fields, which possesses many distinct minima. These 
minima are meta-stable and can decay via bubble nucleation-type tunneling or quantum 
thermal processes. In the early universe, one expects all the minima to be populated as the 
temperature decreases and the fields begin to settle. The minima with the highest values 
of potential energy then begin to infiate, and are expected to rapidly dominate the volume. 
After a few string times the universe will be almost completely full of rapidly infiating regions, 
and as time continues to pass the metastable field values in these regions will begin to decay 
by bubble nucleation. Since we live in a region with very small vacuum energy, we expect 
that our observable universe is contained in such a bubble. 

The most promising way to test this idea is to look at large scale cosmological structure. 
If we model the string theory landscape as a scalar field coupled to gravity with a potential 
with several minima, universes inside a bubble which formed via a tunneling instanton have 
certain characteristic features: they are open, they are curvature dominated at the "big bang" 
(rendering it non- singular), and they have a characteristic spectrum of density perturbations 
generated during the infiationary period that takes place inside them [H O [6l [71 [H HI |10] . 

Moreover, because the bubble forms inside a metastable region, other bubbles will appear 
in the space around it. If these bubbles appear close enough to ours, they will collide with 
it [m US [13]. The two bubbles then "stick" together, separated by a domain wall, and 
the force of the collision releases a pulse of energy which propagates through each of the 
pair. This process (as seen from inside our bubble) makes the universe anisotropic, affects 
the cosmic microwave background and the formation of structure. A collision between two 
bubbles produces a spacetime which can be solved for exactly [TH [13, \W\ . 

In [12], we began the study of this process. We were primarily interested in the effects 
on the cosmic microwave background (CMB). Our solution shows that the CMB sky will be 
divided into two (unequal) halves along a circle. The size of the circle depends on when the 
collision took place. The temperature fluctuations in the disk on the side of the circle away 
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from the collision will be weakly affected (by the gravitational backreaction of the energy 
dumped into our bubble) in a relatively simple way, which we calculated. The disk on the 
other side will be more strongly affected, and in a way that depends in a more complex way 
on the model which we did not calculate. In this paper we will use a simplified model to 
compute these latter effects. 

Our model consists of a collision between our bubble — which we treat as a thin-wall 
Coleman-de Luccia bubble which undergoes a period of inflation after it formed (which we 
model as a de Sitter phase) followed by a period of radiation domination — and another 
different bubble, the details of which are largely unimportant. We assume that the tension 
of the domain wall between the bubbles and the vacuum energy in the other bubble are 
such that the domain wall accelerates away from us. Such collisions can produce observable 
effects which are not in conflict with current data. We will also assume that the infiaton 
is the field that underwent the initial tunneling transition, so that the reheating surface is 
determined by the evolution of that field and can be affected by the collision with the other 
bubble through its non-linear equation of motion, and that no other scalar fields are relevant. 

After some number of efolds of inflation, our bubble will reheat and become radiation 
dominated. The subsequent expansion and the geometry of the reheating surface itself will 
be modified by the presence of the domain wall. The effect we will focus on is the change in 
the reheating surface brought about by the domain wall. 

There are a number of anomalies in current cosmological data, including a cold spot 
[13 [THl Uni Eni El] , hemispherical power asymmetry p2], non-Gaussianities (see e.g. [23]). 
and the "axis of evil" [21]. While the significance of these effects is unclear, there are few 
if any first-principle models which can account for them. Our results indicate that bubble 
collisions naturally create hemispherical power asymmetries spread through the lower CMB 
multipoles and/or produce relatively small cold or hot spots, depending on the time of the 
collision (and hence the size of the disk), and can lead to measurable non-Gaussianities. The 
qualitative features, however, are always the same: a collision with a single other bubble 
produces anisotropics that depend only on one angle — the angle to the vector pointing from 
our location to the center of the other bubble — and which are generically monotonic in the 
angle. While these results are suggestive, further analysis is required to determine which of 
these anomalies can be explained by bubble collisions. 

While we will not study structure formation in this paper, it is interesting to note that 
bubble collisions may also create coherent peculiar velocity flows for large-scale structures 
across large swaths of the sky. This may be interesting in light of the recently published 
results [25] . 

2 The Scenario 

In this section we will set up the basic scenario for the bubble collision we seek to study. The 
details of the dynamics of general collisions were studied in [15] and we refer the reader there 
for a more complete treatment. We will focus on a collision where our bubble is dominated 
by positive vacuum energy when it forms, so that it undergoes a period of inflation, and we 
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will assume the pressures are such that the domain wall is accelerating away from the center 
of our bubble at late times. This second assumption is valid for collisions between two dS 
bubbles if ours has a smaller cosmological constant Q 

Because of the difficulty in finding analytic solutions to the full problem of a scalar 
field with a non-linear potential coupled to gravity, we will make a number of simplifying 
assumptions. We enumerate them here: 

• Unlike in [15], we ignore gravitational backreaction from the shell of radiation released 
into our bubble by the collision and focus instead on the effects on the infiaton's 
evolution and reheating. 

• We approximate the infiaton potential as linear when computing its evolution inside 
our bubble. 

• We allow the field's derivative to be discontinuous across the edge of the lightcone of 
the collision to model the effects of non-linearities. 

• We treat the wall separating our bubble from the false vacuum, and the domain wall 
between our bubble and the other bubble, as Dirichlet (fixed value) boundary surfaces 
for the infiaton. 

• We treat the walls and the edges of the collision lightcone (where the field derivatives 
change) as thin. 

• We assume infiationary perturbations are unaffected by the collision. 

• We model the effects of the collision on the CMB by computing the redshift to the 
reheating surface as a function of angle, rather than evolving the fluctuations in the 
modifled spacetime. 

We will discuss some of these assumptions in more detail later. 
2.1 The Bubble Collision 

In this section we will review the bubble collision spacetime. Readers familiar with the 
details of [15] may wish to skip ahead to Sec. 12.21 

A single Coleman-de Luccia bubble in four dimensions has an 5*0(3, 1) invariance (inher- 
ited from the SO {4) symmetry of the Euclidean instanton [1]). When another bubble collides 
with ours a preferred direction is picked out, breaking the symmetry down to S0{2, 1). As a 
result, one can choose coordinates in such a way that the metric takes the form of a warped 
product between two-dimensional hyperboloids {H2) and a 2D space. 

^It also occurs when our bubble has a larger cosmological constant or in collisions with AdS bubbles when 
certain conditions involving the tension of the wall are satisfied [HI [151 HSl ES] - 
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Under the assumption that all the energy in the spacetime remains confined to thin 
shells or is in the form of vacuum energy, one can find the solution describing the collision 
by patching together general solutions to Einstein's equations in vacuum plus cosmological 
constant with this symmetry p^[T5l[27] . The overall result is that the spacetime inside such 
a bubble is divided into two regions — one that is outside the future lightcone of the collision 
and can be described using an S0{3, 1) invariant metric (a standard open Robertson- Walker 
cosmology), and one which is inside, is affected by backreaction, and is most conveniently 
described using 2D hyperboloids. 



de Sitter space The general 5*0(2, 1) invariant spacetime with a positive cosmological 
constant is given by (see jlSj for a full discussion of the geometry and causal structure) 

ds^ = ^ + g{t)dx^ + fdHl (2.1) 

where 

9it) = 1 + ^ - p (2-2) 

and 

dH^ = dp"^ + sinh^ pdip'^ (2.3) 

is the metric on the unit 2-hyperboloid and x ~ x + ni. We will find for most collisions 
that to '^t at cosmologically relevant times and thus we approximate to ~ 0, giving g{t) ^ 
1 + This is a metric for pure dS space with radius of curvature i. Hereafter, we set 
^ = 1 (it can be easily restored with dimensional analysis). These coordinates (which we 
denote the H2 coordinates) are related to the 5*0(3, 1) ot H3 invariant ones by 



cosh r = Vl + cos 



X, 



sinh r sinh ^ cos 6^ = Vi + t^sinx, 

(2.4) 

sinh T cosh ^ = t cosh p, 
sinh T sinh ^ sin 6 = t sinh p, 

ip = f, 



which gives a metric 



ds"^ = -dr^ + sinh^ rdH^ 
dHl = d^^ + sinh^ ^{de^ + sin'^ edip'^) 



We display the two coordinate systems for dS space in Fig. [H 



Radiation domination We will also need the metric for a radiation dominated (RD) 
universe. The general 5*0(3, 1) invariant Friedmann-Robertson- Walker (FRW) metric is 
given by 

ds^ = -dr^ + a^{T)dHl. (2.6) 
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Figure 1: Plot showing the two coordinates systems for dS space. The straight (red) hnes are 
t,x, and the curvy (blue) lines are r, ^. The values on the axes are in conformal coordinates 
where the vertical axis, tan~^ t, runs from to 7r/2 and the horizontal axis, x, from — 7r/2 to 
7r/2 (with£= 1). 

The transformation to H2 coordinates is 

a(r) sinh^cos6' = z, 

a(r) sinhf sin 6' = tsinhp, 

.X (2-7) 
a(r) cosh ^ = tcoshp. 

For radiation domination we have a(r) = V Cr, with C a constant. This gives a metric in 
H2 coordinates 



-4t^ + + 4t2) ^2 ^2^2 _ 4^2 (^2 _ 



'^^^ [A{e-z^)-c^]dtdz + edHi 



(2.8) 



C2(t2 _ ^2 



The radiation dominated part of the spacetime will be affected by the domain wall and the 
collision. Because of the reduced symmetry (relative to de Sitter), we are currently unable 
to write down the most general radiation dominated solution with 5*0(2, 1) invarianceH 
However, we expect that at least in a broad class of scenarios the effects on the CMB due to 
backreaction will be subleading to the other effects we will study, and so we will use this form 
for now. It would be interesting to consider these metric effects in the future, for example 
by solving the equations numerically. Such a study is currently in progress [28] . 



2.2 The Reheating Surface 

With these building blocks we can construct our model. Our universe appears as a de Sitter 
bubble that collides with another bubble, forming a domain wall along the interface. After 

^The problem is equivalent to finding the smooth metric describing a Schwarzschild black hole embedded 
in an asymptotically FRW spacetime. 
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Observer Radiation Shock 




Bubble Walls 



Figure 2: Sketch of the colhsion scenario. Due to the domain wall, the shape of the reheating 
surface deviates from that of the no collision case (displayed in the long dashed line). 

the initial period of inflation the universe reheats and transitions to radiation domination. 
We model this transition as sharp, matching the metrics and Hubble constant across the 
reheating surface. We ignore matter domination in calculating the effects of the collision on 
the CMB, and further assume that the effects on the surface of last scattering are related 
to the effects on the reheating surface simply by the redshift along "free-streaming" null 
geodesies. In other words, we assume the angular dependence of the redshift to the last 
scattering surface due to the effects of the collision is given by the angular dependence of the 
redshift to the modified reheating surface. Ignoring the evolution of perturbations between 
reheating and last scattering is a valid approximation on angular scales larger than that 
of the first acoustic peak in the CMB (about one degree). On smaller scales there will be 
a "smearing" due to both this evolution and to the finite thickness of the last scattering 
surface. Due to these effects, we will focus on disks considerably larger than this size. We 
will also ignore the effects of the collision on the perturbations generated during inflation. 
While this should be justified when the overall effects are small enough to be consistent with 
observation, it would be interesting to study this more carefully. To summarize, in our setup, 
the main effect is to modulate the standard inflationary perturbations at last scattering by 
a function dependent on the angle to the collision direction. We sketch the scenario in Fig. 

m 

We work in the thin wall approximation throughout, treating the transition between the 
two parts of the sky (corresponding to the regions of last scattering which are inside and 
outside the collision lightcone) as sharp. It would be interesting to extend the present work 
to thick wall systems as effects such as the finite thickness of the domain walls and the shell 
of radiation emitted from the collision point could play interesting roles. 

In a standard Robert son- Walker open universe, reheating occurs on a surface of almost 
constant density that is 5*0(3,1) invariant (i.e., hyperbolic 3-space). During inflation the 
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(a) (b) 

Figure 3: The scalar field transitions from some vacuum to ours along both the domain and 
bubble walls. It is possible for the field to transition from the same direction along both (a), 
or from different directions for each (b). 

surfaces of constant energy density are also surfaces of constant value for the inflaton field, 
and the reheating surface corresponds to the point in the inflaton potential where slow roll 
comes to an end. This surface defines a preferred reference frame, since the radiation and 
matter fluids will be at rest on average with respect to it. In our scenario, the collision 
and domain wall break this symmetry. There is no longer a preferred frame and thus the 
standard notion of "comoving" no longer holds. There remains a remnant of this symmetry: 
a two-dimensional set of comoving observers related to each other by transformations which 
leave the two-dimensional hyperboloids invariant. Observers outside this plane are no longer 
equivalent. 

We define the reheating surface in our scenario as a surface of constant scalar field, 
4>{x,t) = 00- The presence of the domain wall will alter the shape of this surface in a 
non-trivial way, breaking isotropy. To fully solve this problem we would need to solve the 
non-linear coupled scalar field and gravity system for the collision spacetime. Instead we 
will make an approximation. We expect the scalar field to have a set of values 0^ in the 
other bubble, some other set of values (pn in our bubble and a third set of values in the 
metastable vacuum outside both bubbles (pM- In the thin wall limit the scalar field jumps 
from near some other local minimum (either the false vacuum our bubble tunneled from, or 
the vacuum inside the bubble it collided with) to a location a bit up the "hill" from our local 
minimum. This transition occurs all along the domain and bubble walls, so the scalar field 
is constant along each of these surfaces. 

Without loss of generality we choose = along the bubble walls. At the domain wall, 
we have two options. The field could transition at the collision and domain wall from the 
same direction in field space as the tunneling that initially created our bubble Fig. [3](a). In 
this case, to a good approximation we have = along the domain wall as well. However, 
it is also possible for the field to transition from the other side of the potential Fig. [3](b), in 
which case we have (p = k ^ along the domain walljfl 

■^In the string theory landscape we expect the scalar field space to be multi-dimensional, in which case 
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The simplest potential for cf) which allows it to roll inside our bubble is linear: V{(f)) = fi(j), 
which can be thought of as the first term in an expansion around the exact potential on the 
inflationary slope. Our first task is to solve for the surfaces of constant in this geometry 
and with these boundary conditions. 



2.3 Observer Horizons 

Before proceeding, it is instructive to determine how much of the reheating surface is inside 
the past lightcone of an observer at some given time. We neglect the period of matter 
domination, which should not strongly affect the conclusions. We wish to calculate the 
radius C,r of the surface defined by the intersection of the observer's past light cone with 
reheating, from which we can determine the corresponding values of the coordinates z and x 
(see (12. 4p and (12. 7p ). The past lightcone of an observer at the origin ^ = consists of a set of 
geodesies which in coordinates remain at constant 6 and (p. In the radiation dominated 
spacetime the null geodesic is given by 



dr 



(2.9) 



Ae = -^A(v^)- (2.10) 



We want to match this up to the dS space. At reheating we need to match both the scale 
factor a and Hubble constant H of the two metrics. This gives 



a = ^/ CTr = sinhxr, (2-11) 
a 1 1 



a 2Tr. tanh 



(2.12) 



where Tr is the time of reheating in the radiation dominated spacetime, and fv is the 
same for the dS spacetime (we could make the coordinate r continuous at the reheating 
with a shift, but it will be more convenient not to). We will use tildes for the dS time 
throughout this subsection; the coordinates on the itself are continuous. The solutions 
are 

1 

Tr = -tanhr^, 

C = 2 sinh^r cosh f^. (2-13) 
Using the fact that temperature redshifts as 1/a, 

e^* = Tr/Tn = ttn/ttr = a/t"/^, (2.14) 

where n subscripts refer to the observer's time and this defines A^,=. 
Hence, r„ = e^^'Tr. This gives 

= ~ 2e^*-^'- = 2e^*-^, (2.15) 

cosh 



the field could come in from the "side" . It would be interesting to investigate how this affects our results. 
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where N is the number of efolds of inflation and we have used A^, A^=k ^ 1. In this case 
$,r <^ 1 and the flatness problem is solved (because ^ ~ 1 is one curvature radius of the H^). 
Hence A^* is roughly the number of efolds required to solve the flatness problem. A simple 
calculation demonstrates that this same condition is required to solve the horizon problem. 

In [15| we calculated the location of the intersection of the future lightcone of the collision 
with the reheating surface, and found that the parameter can be roughly as large as t^, 
where tc is the time of collision and to is the backreaction parameter appearing in (12. ip . 
We found that t^. could be exponentially large and still have the collision within the ^ = 
observer's past lightcone. This was done in a toy model by patching a curvature dominated 
open universe onto the dS space at reheating. One can see from the calculation above that 
patching on radiation domination rather than flat space dramatically changes this estimate, 
and that tc cannot be exponentially large and still be visible to an observer at ^ = 0. This 
justifies taking ~ 0, as we have done here. 

Because these collisions break part of the isometry group of open Robertson- Walker 
cosmology, they define a preferred "center of mass" frame in which the two bubbles nucleate 
simultaneously. In our coordinates, observers at rest in this frame move on geodesies of 
constant z (or x). The intersection of this set of observers with the set comoving inside 
the bubble cosmology are geodesies with z = ^ = 0, spanned by the H2. In the analysis 
above we focused on such a "center of mass" (COM) observer. Observers who are not at rest 
in this frame (which we will refer to as "boosted observers"), including those who formed 
from matter comoving with respect to the unperturbed part of the reheating surface, will see 
effects which depend on their boost relative to the COM observers. Those boosted towards 
the collision can see bubbles that nucleated farther away in the COM frame (and hence 
collided with ours later). For a given scenario, observers boosted towards the collision will 
see effects enhanced relative to the COM observer (essentially because of Doppler shift). The 
opposite applies to observers boosted away. 

We will discuss this in more detail in Sections 14.21 and O 
3 Finding The Reheating Surface 

In this section we'll solve for the geometry of the reheating surface. The computation involves 
several steps. First, we will outline the geometry, causal structure and domain wall trajectory 
in our bubble. These will determine the boundary conditions on the scalar field. We then 
solve the scalar field equation for a linear potential with boundary conditions that it is 
constant on the bubble and domain walls. Armed with the solution for we can find the 
reheating surface by solving (p{z,t) = (po, where 0o is the field value where the field ceases 
to slowly roll and reheating begins. 
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3.1 The Domain Wall Trajectory 



The trajectory of the domain wall in x,t (recall that the wall is homogeneous on the H2) 
was found in 115] 



.2 W ^ 



x' = —- [-1 + —-] , (3.1) 

git) V git) J ^ ' 

where here a dot denotes differentiation with respect to the proper time s along the wall. At 
late times the solution for t{s) is 

tis) = Toe"^ (3.2) 

To a constant. If we assume the other bubble colliding with us is also a dS bubble with 
radius of curvature £2, then 



A^ = l + -.U^ + -^-l (3.3) 



(in units of the radius of curvature of our bubble). At late times we can use g{t) ^ and 
plugging in the solution for t{s) we fincfl 



-W = ,/l-^(i-i)+:^.-«Q-y+:r„ (3.4) 

where (tc, Xc) are the coordinates of the collision. < a < 1 is related to the proper 
acceleration of the wall, with a = corresponding to no acceleration and a = 1 to infinite 
(so that the trajectory with a = 1 is null). 

3.2 Null Rays And The Structure Of H2 dS Space 

In this subsection, we'll set up the geometry and boundary conditions we'll need to solve the 
scalar field equation in the collision background. For dS space, the null rays are given by 

doc 1 

— = ±- ^ x(t) = ±tan~4 + A;, (3.5) 
dt g 

where k is an integration constant. There are three null surfaces important to us: the 
two walls of the initial bubble and the lightcone of the collision. We choose to orient our 
coordinates such that our initial bubble is nucleated at {t, x) = (0, 0). Running through each 
line we find 

right wall — x = tan~^ t, 
left wall — > x = — tan~^t, (3.6) 
radiation line — > x = tan^^ t — 2 tan~^ tc, 

and the collision point is at {tc,Xc) = (tc? — tan^^t^)- 



*It is straightforward to generalize this to a coUision where the other bubble is AdS, results will be similar. 
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x=-tanH l/t-7i/2 



Figure 4: Sketch of the coUision scenario with the null lines labeled and approximate solutions 
for late times. 

We'll also need these rays for large t. We can take the asymptotic expansions of (13. 6p to 
gelEl 

right wall — > X = 1 — , 

1 TT 

left wall — >■ ^ ~ ~ 2' ' i'^-'^) 

V . 1 2 vr 

radiation line — »• x = 1 , 

t 2' 

and the collision point is at {tc,Xc) = {tc,j- — f)- The geometry and causal structure are 
summarized in Fig. HI Note that region I is outside the lightcone of the collision, and the 
geometry can effectively be described in coordinates up to the radiation line. 

3.3 General Solutions For The Scalar Field 

In this subsection we'll find the general solutions to the scalar field equation (with a linear 
potential) in the regions before and after the collision. Let's begin by looking at region I. 
Here we can make use of the coordinates. We have 

□0 = ^-^dr (sinh^ Tdr(t){T)) = fi, (3.8) 

sinh T 

where /i is the coefficient of the linear term in the potential. We require that (j) and its 
derivative vanish along the bubble walls at r = 0. The general solution satisfying these 

^We have here also expanded for large t^, for smaller values of tc we can convert back by exchanging 
1/tc ^ — tan"-'^ tc + tt/2. 
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boundary conditions is 



3 



ln-^^--tanh2 (r/2) 



(3.9) 



At large and small r this becomes 



0(r > 1) ^ -- Insinhr ^ -- In(tcosx) (3.10) 
3 3 

2 

0(r<l) ^ ^ (^Vl + t2 cosx - l) , (3.11) 

where we have used the coordinate transformations (12.40 . 

In region II we solve for using the H2 coordinates, as it is in these coordinates that the 
boundary conditions on the domain wall and collision lightcone are simple. Since the field 
is constant on the H2 we have 

□0 = ~dt + e)dtct>] + Y^dl^ = (3.12) 

The general solution is 

(j){t,x) = /(a; - tan"4) - ^/'(x - tan~U) +5((a; + tan~4) + ^5('(a; + tan"4) 

- ^ln(l + t^), (3.13) 

where / and g are arbitrary functions and the primes denote derivatives with respect to the 
argument of the respective function. At large t, this is given by 

0(t, x) ^ fix + 1/t) - J fix + 1/t) + g{x - 1/t) + ^g'{x - 1/t) - | In t. (3.14) 



3.4 Matching Solutions 



We have the general solution for the scalar field in the H2 symmetric background, as well 
as the solution in region I fl3.10p . before the collision. To find the particular solution for 
(f) in region II we make use of two boundary conditions: the field must be continuous at 
the radiation line, and it must equal a constant at the domain wall, (pidomain = k. Since 
reheating occurs after inflation we can use the late time approximations to the general 
solutions. Carrying out these steps is somewhat technical; the details are in appendix |Al 

In Fig. [5] we draw an example with A^^, = 60, tc = l,a = 0.3, /i = 1, and k = 0. The 
presence of the domain wall causes the surfaces to "curl up" , eventually going null and then 
timelike. That part of the spacetime is very strongly perturbed by the collision — inflation 
and reheating will not occur there, and the linear techniques used in this paper do not apply. 
Since these regions are highly perturbed away from a standard cosmology, the effects would 
be in drastic conflict with current observations, and thus we will assume that this region is 
well outside our horizon. 
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Figure 5: Plot of surfaces of constant (f) at late times for A^* = 60, a = 0.3, tc = 1, = 1, 
and k = 0. The blue lines are the surfaces if no collision had occurred, while the red lines 
are the actual surfaces. Black lines are the bubble walls (projected forward as if there were 
no collision) and the green line is the domain wall. This plot is in the same coordinates as 

Fig. m 




Figure 6: Plot of surfaces of constant (f) at late times for A^^, = 60, a = 0.3, k = 0, jj = 1, 
and tc = 3. The blue lines are the surfaces if no collision had occurred, while the red lines 
are the actual surfaces. Black lines are the bubble walls (projected forward as if there were 
no collision) and the green line is the domain wall. 
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Figure 7: Plot of surfaces of constant (p at late times for A^,,, = 60, a = 0.5, tc = 1, = 1, 
and k = —1. The blue lines are the surfaces if no collision had occurred, while the red lines 
are the actual surfaces. Black lines are the bubble walls (projected forward as if there were 
no collision) and the green line is the domain wall. 

We can make some other general statements. Increasing the collision time tc pushes the 
difference between the actual surface and that of no collision farther out in x. In Fig. E] we 
draw the same scenario as in Fig. O except with tc = 3. The smaller a is, the more rapidly 
the surfaces curl up. Taking a — 1 in the case A; = removes the effects of the collision 
entirely (which is a good check of our numerics). 

For A; < we find that the reheating surface bends down relative to the surface of the 
no collision scenario. In Fig. [7] we draw a case with A^^, = 60, tc = l,a = 0.5, k = —1, and 
fi = 1. Increasing tc causes less of the sky to be covered by the collision, while the smaller a 
gets the more pronounced and closer to the observer the effects are. Taking A; > is similar 
to the scenarios with k = except that the repulsion from the domain wall will be increased. 



In this section we will use our solution for the modified reheating surface to find the redshift 
from it to observers in the bubble. Using the results of Sec. 12.31 we find the projection 
of the observer's past lightcone on the reheating surface and calculate the redshift as a 
function of angle, taking into account the doppler shift due to the fact that the reheating 
surface is not comoving. We will again begin by considering COM {z = 0) observers, 
and later generalize to boosted observers. In doing this analysis we will ignore the effects 
of gravitational backreaction from the collision and modified reheating surface, as well as 
the modification to a(r) from spatial curvature . These approximations are justified when 
^ 1, at least for observers away from the part of the spacetime that is strongly modified. 



4 The Redshift 
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4.1 Null Geodesies In Radiation Dominated Spacetime 



Starting from the results of Sec. 12. 3[ we can express the geodesies in H2 coordinates using 
the coordinate transformations (12. 7p 

t = (^^{^r - + V^) ^cosh^e-sinh^esin^^, 

z = (^^(^,-0 + v^^sinh^cos^, (4.1) 

tanh p = tanh ^ sin 6, 

where for each geodesic 6' is a constant and the set of all geodesies is given by 6 taking all 
values from to tt, with 

C ^ -e^^, Tr ^ -. (4.2) 
2 2 ^ ^ 

Since Zn = the observer's time t„ is given by 



1 ^.r ^.r 1 



tn = VC^a ^ \l -e2^e2^* ^ -6^+^* . (4.3) 



Using this we can solve for by setting = 0) = t„, giving 



= 4e-2^ ( _ le^ ) ^ 2e^*-^. (4.4) 



2 



Putting this all together we get 



t ^ \ (e^*+^ + - ^e^^ j ^cosh^ i - sinh^ i sin^ ^, 

^ ^(^e^*+^ + e^-|e2^^sinhecos^. (4.5) 

To compute the redshift we have to take into account the relative motion of the reheating 
surface and observer. To do so we need another piece of information about the geodesies, 
namely their 4-momenta at both the time of emission and when they intersect the observer. 
We can again make use of the i^a coordinates. Differentiating both sides of the z equation 
of (12. 7p with respect to the affine parameter we find 



z = jf 



' Ct cos ^ ^^^^ ^ + cosh ^ , (4.6) 



where p" = r and p^ = ^. Since these geodesies have constant 9, p^ = g^^p^ is conserved. 
From (12. 9p . p'^ = —y/Crp^. We can solve for in a similar way to get 

J. Cz tcoshp 

P = - 2(t^-.2)3/2 PC + 7^PC^Osg, 

f Ct 2; cosh p y ■ ) 
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4.2 Boosted Observers 



So far we have considered the special case of observers at rest in the frame in which the 
bubbles nucleated simultaneously. In our coordinates these COM observers are those in the 
z — plane, and can always be taken to be at the origin of the H2 {p = 0) by a shift in the 
origin of coordinates. 

An observer comoving in the open universe inside the bubble is not in general at rest 
in this frame. But since the boosts "parallel" to the domain wall are trivial {i.e. one can 

always choose coordinates so any comoving geodesic lies in the t, z plane and at the origin 
of the H2), the most general comoving observer is related to an observer at ^ = 0, due to 
the ifs isometry, by a "boost" in the z direction only. The isometry is 

r' = T (4.8) 

cosh ^' = 7 (cosh ^ — 13 sinh ^ cos 9) 
sinh ^' cos ^' = 7 (sinh ^ cos 6 — j3 cosh ^ ) , 

where 7 = (1 — /?^)~^/^ is a "boost" around the nucleation point of our bubble. (3 > Q 
corresponds to an observer closer to the collision (boosted to the left in our diagrams). In 
H2 coordinates this becomes 

t'coshp = 'J {tcoshp - l3z) (4.9) 
z' — J {z — pt cosh p) 

It is clear that a boosted observer will not see the same redshift pattern from the collision 
in the CMB as a COM observer, even at equal comoving time. To compute this effect we need 
to repeat the calculation we have just done, only now for an observer following the boosted 
trajectory. We can find the boosted lightcone by acting on the lightcone of the z = observer 
with the boost (this works since the boost is an isometry of the unperturbed spacetime). 
For brevity we will not include the details here; the calculation is straightforward. 

One should be careful to distinguish two cases: structures whose comoving trajectories 
intersect the reheating surface in the region outside the hghtcone of the coUision and only 
pass into it much later, versus those that do not. We expect the first class of objects 
to be comoving in the unperturbed spacetime (up to peculiar velocities), but the second 
class formed in a part of the spacetime where the reheating surface was perturbed. These 
structures may be locally at rest with respect to the frame defined by the modified reheating 
surface, and will possess a coherent "peculiar velocity" with respect to the first class. Our 
analysis here is valid only for observers in the first class, but in the regions of the space 
where the effects are small it should be approximately correct. 

4.3 Finding The Redshift 

Armed with the geodesies and the reheating surface, we have all the information we need 
to compute the redshift as a function of the angle on an observer's sky. For a given angle 
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Figure 8: Relative redshift (rs/ro) versus cos^ for the scenario (tc = l,a = 0.5, N — N^, 
0,/i = 1). Relative redshift greater than one indicates higher relative temperature. 



5,k 



6, if on the observers sky, we follow the null geodesic back to where it intersects the modified 
reheating surface. We take the dot product of the geodesic's 4-momentum with the unit 
normal to the reheating surface at the intersection point, and then the dot product with 
the tangent vector to the observer's geodesic. The ratio of these two quantities gives us the 
redshift in the approximation we mentioned before. 

We first calculate the reheating surface and normals from our solution for (f) in Sec. Ei 
The surface is found by solving (p{t,x) = 0o for t as a function of x and 0o- For a given 
potential V{(f)), 0o determines the number of inflationary e-folds N. In the end we have the 
modified reheating surface for a given set of parameters (A^, iV*, tc, a, /i, k). 

Once we have the intersection point, we can compute the vectors normal to the surface 
by taking the exterior derivative of (f){t, z) (in the RD part of the spacetime). We normalize 
these using the RD metric fl2.8p . The tangent 1-form along the observer's trajectory is 
proportional to dt for an unboosted observer at ^ = 0, and to dr ~ tdt — zdz in general. The 
redshift is given by 



En _ -{P^n^)n 



Er 



— {p^n 



(4.10) 



II Jr 



where n refers to the observer's position, r to the point on the reheating surface, and are 
the unit 1-forms at each. Recall if there was no collision then the redshift would be 



To 



(4.11) 



which is also the redshift for any geodesies that originate from parts of the reheating surface 
that lie outside the future lightcone of the collision. In Fig. [8] we show some results for the 
scenario (tc = l,a = 0.5, N — N^, = 5, = 0,/i = 1). We plot r^/ro as a function of cos^ 
where < < vr is the angle on the sky with 6 = it pointing towards the collision (the 
redshift s are axially symmetric). 

We can make some general observations about the redshift. The modified reheating 
surface leads to bluer photons when k > 0, the bluest coming from the angle closest to the 
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Figure 9: Relative redshift (rs/ro) versus cos^ for the scenario (tc = 1.05, a = 0.1, N — = 

3,k = —1, = 1). 

domain wall. This makes sense, as the reheating surface is closer to the observer than it 
would be without the collision. However, the doppler shift from the motion of the reheating 
surface with respect to the observer also plays a role. The effect is primarily a dipole, but 
only on the range of angles for which null geodesies intersect the modified reheating surface. 
Once we are viewing angles which correspond to parts of the reheating surface outside the 
future lightcone of the collision, the redshift is tq. The smaller a is, the more pronounced the 
effects of the collision are, as the domain wall is closer to the observer and accelerating away 
more slowly. Raising tc means less of the sky is within the future lightcone of the collision. 
Thus, there are a wide range of values for both the maximum blueshift and the angular size 
of the part of the sky affected by the dipole piece. If tc is much larger than 1 in Hubble units 
the COM observer is likely not within the future lightcone of the collision (as we mentioned 
in Sec. [23D . 

If we look at scenarios where k < we see that the modified reheating surface will lead 
to slightly redder photons, again the reddest coming from the angle closest to the collision 
direction. This is due to the reheating occurring earlier in that direction. In Fig. we 
display a typical scenario with [tc = 1.05, a = 0.1, N — N^, = 3, k = —1, fi = 1). 

Let's see how this effect varies with the number of e-folds N. In Fig. [TD]we show some 
sample values for a collision with tc = 1, a = 0.5, = 0,/i = 1 at angle cos^ = —0.9. We 
plot it for 1 < — A'^* < 11 and draw a best fit line. The best fit line at this angle is 
r^/ro ~ 1 + 0.77 e^*~^, which decreases exponentially in N as expected from inflation. We 
expect that the general behavior will go as 

-{a, tc, 9,k,fx)^l + tc, 9, k, /i)e^*-^, (4.12) 

where < |^| < 1 is a function of the specifics of the collision and the angle. Thus, unboosted 
observers can detect effects from the collision so long as the number of e-folds of inflation is 
not too large. 

Observers boosted towards the collision will see more of the modified part of the space, 
and hence the angular size of the affected spot will be larger. On the other hand, observers 
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Figure 10: Relative redshift {vs/ro) versus N — N^, for the scenario (tc = l,a = 0.5, k = 
0, n = 1) at angle cos^ = —0.9. The dots are actual values and the line is a best fit line with 
equation r^/ro = 1 + 0.77 e^*~^. 

boosted away see the opposite (smaller, weaker spots). For both, the angular dependence 
within the spot is similar in form to the COM observer, constant away from the collision 
and linear towards the collision. It is interesting to note that even long inflation cannot 
completely remove the effects on the skies of some boosted observers, since there are always 
parts of the space near the domain wall. However these observers are not well-approximated 
by the treatment carried out here, since in the limit of large boost and long inflation they 
formed from strongly affected parts of the reheating surface. Nevertheless, it is possible that 
such observers see a more or less standard cosmology. The study of this question is beyond 
the reach of this paper. 

tc larger than of order 1 will take the COM observer outside the future lightcone of the 
collision. However, it is again true that a sufficiently large boost towards the collision will 
bring the observer back into the collision's lightcone. Hence, no matter what the values of tc 
and are, some set of observers will always be able to see the effects of the collision. How 
probable we are to find ourselves in that set is again beyond the scope of this paper, but see 

m- 

5 A Simple Model 

The preceding results show that the most natural way to describe the signature of a collision 
is in position space, with a simple piecewise linear behavior for the redshift as a function of 
cos 6. However, in order to make direct comparisons to previous CMB analyses, we will use 
this result to explore the consequences for the CMB multipoles. Looking at Figures [8] and [9] 
we see that the redshift can be approximately modeled as 



where tq is a constant that gives the redshift with no collision, and x = cos 6. The param- 
eter M gives the maximum value of the relative blue or redshift due to the collision. The 




(5.1) 
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Figure 11: An example of the temperature (excluding fluctuations) as a function of angle in 
our toy model. 



parameter xt is the cosine of the angle where the forward light cone of the collision ends, 
i.e. for X > Xt we are looking at the unperturbed reheating surface. 

Including the density fluctuations at reheating, the temperature today is 

r(n)=T^r(n)(l + 5(n)), (5.2) 

where Tq is a constant that sets the overall temperature and (5(n) are the temperature 
anisotropics from inflationary density perturbations. We will take 5(n) to be Gaussian 
random variables with zero mean. We plot an example of the function TQr(n) in Fig. [TTl 

To flnd the CMB multipoles we need the fluctuations of the temperature about its average 
value. This is given by 



5T{n) ^ roV(n)(l + ,5(n))-To 
To To 

where 



(5.3) 



(5.4) 



is the average temperature and we've used the fact that J (Pn 5{n) = 0. 

Several authors have considered anisotropic CMB spectra where the temperature fluctu- 
ations are modulated by a function [SOj [22l [31]. We can parallel some of the formalism in 
[30] with some modiflcations. We write 

^ = A{n) + /(I + t.(n))i?(n), (5.5) 

where B{n) = 6{n), f = T^/Tq and 

w{n) = 6(xr — x) {xt — x) ^ . (5.6) 
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In [30], A{h) is a random variable but with a distribution that differs from that of 5(n). 
In our case, we want to find the expected CMB multipoles with fixed parameters for the 
colhsion; i.e., we will average over density fluctuations while holding fixed the collision 
redshift function r(n). Ideally, we would use a microscopic model (for example, the string 
theory landscape) to compute the probability distribution for the parameters tc and (3 — 
which are determined by the location and time of nucleation of the other bubble and are 
random, since the nucleation is a quantum event — and hence a probability distribution for 
r(n). However, without detailed knowledge of the probability measure for collisions we 
are not currently equipped to carry out this procedure. In fact, there could be multiple 
collisions in the past lightcone of the observer, each with an associated time and position for 
its nucleation. It would interesting to study this, but here we will proceed to study a single 
collision with fixed values for the parameters. 

A{n) is given by 

^(ft) ^ 5^^^^ - 1. (5.7) 

We can expand the fluctuations in multipoles via spherical harmonics. 

6T{h) 



Tn 



Lm 



For details on our conventions for the spherical harmonics see appendix [Bl We can expand 
A, B and w in multipoles as well 

A{ri) = y^^aimYim{n), 

lm 

B{n) = ^6,^ri^(n), (5.9) 

lm 

W{h) = ^Wi„Xlmir^)- 



lm 

Since we have chosen coordinates such that 6 = n corresponds to the collision direction, 
both w and A have only m = multipoles. Other coordinates can be obtained via a 
straightforward rotation using Euler angles. 

We are now in a position to write the multipoles of the temperature fluctuations. The 
coupling between w(n) and B{h) becomes a convolution in fourier space 

tlm = aioSmO + fkm + R\mKm, (5.10) 

hh 

where the coupling matrix is written in terms of Wigner 3j symbols 



oW._. (2/ + l)(2/i + l)(2/, + l) //, h l\(h k ^ ....^ 



where wi = wiq. 
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The functions A[n) and w{h) have non-zero 1-point functions: 




(5.12) 



while 



(blm) = 0. 



(5.13) 



Thus we can quickly compute the multipoles for the one-point function 




(5.14) 



The only anisotropy we see at this level is a result of the collision itself. Consistency with 
the observed CMB dipole requires that \M — 1| < 10~^. 

A more useful object is the 2-point function. We have 



Ci is the two point function in the absence of a collision. We'll discuss this in a bit. 
Squaring (15.101) and taking ensemble averages using the above rules we get 



This determines the two-point function in terms of and the collision redshift function. 
As a first approximation we can assume that comes from primordial fluctuations that 
are unaffected by the collision. We use a spectrum generated with CMBFAST code [32[ |33] 
using concordance cosmology values from WMAP [Ml [35] . 

5.1 Results 

In this section we will present our results for a few simulations of our toy model. The effects 
are predominantly on large angular scales, so we focus on the lowest set of / modes. We will 
look at two examples, one where the collision takes up most of the sky, and another where it 
covers only a small fraction. For each scenario we present two figures. One will display the 
power in each / mode compared to a scenario where no collision occurs, the other the power 
asymmetry as a function of / between the two hemispheres defined by the collision direction. 

The Ci's are 




(5.15) 




m=l 



(5.17) 
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Figure 12: C//Cj°^ vs. / for the scenario M - 1 = 2 x 10-^ xt = 0.8 
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Figure 13: vs. / for the scenario M - 1 = 2 x 10-^ xt = 0.8 



Due to the fact that the aim's only contribute to the m = mode of each multipole, these 
terms gain more power from the colhsion than the modes with m ^ do. This is due to the 
fact that the effects of the colhsion are azimuthally symmetric. 

We compute the power asymmetry between the left and right hemispheres using a crude 
version of the Gabor transform method [36] with a top hat window function centered in 
either hemisphere. There will be effects due to ringing from the edges of the top hat, but 
these don't greatly affect the lower multipoles. 

For our first example, we choose parameters M — 1 = 2 x 10^^ and xt = 0.8. The angular 
radius of the collision is 143 degrees. Note that the total temperature difference between the 
two poles is well below the dipole due to the Earth's peculiar motion. In Fig. [12] we plot 



C;/c/°^ for this scenario. As expected, since the collision takes up a large portion of the sky. 



the lowest multipoles are the most affected, 
asymmetry. 



In Fig. [12] we show the hemispherical power 



In our second scenario the collision takes up only a small part of the sky, and the collision 
causes a redshift within the disk. We choose M — 1 = —7.3 x 10~^ and xt = —0.99, giving 
an angular radius of about 6t = 180° — cos^^ xt = 7° . With these parameters our sky has a 
cold spot, with the coldest point about 20/iK cooler than the average and the effect falling off 
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Figure 14: Ci/cf^ vs. / for the scenario M - 1 = -7.3 x 10-^ xt = -.99. 
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Figure 15: c\^f^/Cl'^^* vs. / for the scenario M - 1 = -7.3 x 10"^ xt = -.99. 

with radius from the center. We have chosen the angular radius and the temperature profile 
to approximate the WMAP cold spot smoothed over a scale of ~ 5° p7l [T8l [T9l[2nl[2T] . 

In Fig. [13] one can see that the power is shifted slightly towards higher multipoles, with 
a first peak in excess power around I = 17. There will be a smaller secondary peak at the 
first harmonic. The quadrapole receives only a small boost in power, reducing it relative to 
nearby multipoles. 

In Fig. [15] we show the hemispherical power asymmetry, which is nearly identical to 
the previous figure (since the effects of the collision are now localized wholly in the left 
hemisphere). We see that collisions of this type can lead to measurable hemispherical power 
asymmetries. 

We have presented results for two specific choices of parameters. The results for general 
parameters all have certain features in common: most of the excess power is in the m = 
modes, the temperature function is monotonic with angle within the affected disk, the power 
spectra are affected primarily in the I modes corresponding to the size of the disk and its 
harmonics, and there is a hemispherical power asymmetry with a magnitude that depends 
on the size and intensity of the disk. 
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Non-Gaussianities: The effects of the bubble colhsion induce sizable non-Gaussianities 
in the CMB temperature map. However, at least for the two examples presented here the 
non-Gaussianity is significant only at low /. In particular, is very small, and y^^''^*^''^' jg 
strongly dependent on /: it oscillates in sign, reaches a maximum near the I corresponding to 
the disk size, and damps rapidly (with an envelope similar to the plots of the excess power in 
the 2-point function). We suspect that an analysis of the WMAP data would not constrain 
these models significantly. We will mention a better statistical test in the conclusions. 



6 Conclusions 

In this paper we have computed the effects of a cosmological bubble collision on the cosmic 
microwave background, using a variety of approximations. Previously we focused on a range 
of parameters where gravitational backreaction significantly affects the CMB even outside 
the collision disk [15]. However this effect is probably too weak to be observable in a more 
reahstic model. Instead, in this paper we concentrate on what we expect to be a stronger 
effect which arises from the fact that the reheating surface — defined by a constant value of 
the tunneling scalar field — is perturbed from the case with no collision. Because of this the 
redshift back to the last scattering surface depends in a distinct and specific way on the angle 
from the collision direction (but not on the azimuthal angle). In particular, our calculations 
demonstrate that this has the behavior seen in Fig. [H] — no effect for cos 6 > cos(7r — 6t) 
and linear in cos 6 for cos 6 < cos(7r — 6t) (where 6t is the angular radius of the affected disk 
and depends on the collision parameters). This form of anisotropy is seen for all observers 
in the part of the spacetime where our approximations are valid, but different observers will 
measure different values for 6t and the temperature gradient inside the disk. 

To compare to the CMB maps, we assume that the inflationary perturbations are unaf- 
fected by the collision and calculate the temperature perturbation map as well as the two 
point function for fluctuations. The deviations in the C;'s are primarily at large scales (low 
/), and depending on the transition angle, affect Z's at harmonics of tt/Ot- Non-Gaussianities 
can be relatively large depending on the parameters, and are best detected using equilateral 
triangles. 

It is possible that this scenario can explain some of the current cosmological anomalies, 
but the analysis we have done so far is not detailed enough to provide a "smoking gun". 
A good statistic for the effect described in this paper would be to pick an origin, measure 
the average CMB temperature in a disk of some angular size around that origin, and then 
find the point and angular radius which maximizes the effect. We have not performed this 
analysis, but the results would be very interesting. 

Many open questions remain. It should be possible to calculate how the inflationary 
perturbations are affected by the collision and include this effect. Aside from the CMB, 
there should be anisotropics in other cosmological observables, such as large scale structure. 
As we mentioned in the introduction, one expects coherent peculiar velocities due to the 
fact that the affected part of the reheating surface is not comoving in the frame of the 
unaffected region, as observed in Other potential signals could be detected using 21 
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cm observations, CMB polarization [37], and better precision on measurements of spatial 
curvature. Finally, it is important to estimate how many bubble collisions to expect in our 
past lightcone, and what the size distribution of disks should be on the sky. 
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A Solving For The Scalar Field 

In this appendix we will describe how to solve for the scalar field in the collision geometry. 
In region I we already have the solution (13.90 . In region II we have to match our solution to 
the region I solution at the radiation line 

Mt, X = 2/t, - l/t - 7r/2) = /(2/t, - 7r/2) - -^f{2/U - vr/2) + g{-2/t + 2/t, - 7r/2) 

t 6 [A.V 

= --Int - ^ln[cos(-l/t + 2/tc - 7r/2)] 
3 3 

= (f)i{t,x = 2/tc-l/t-n/2). 
This is a differential equation for g. We can write u = —2/t + 2/tc — 7r/2 and get 

= -|ln[cos(V2 + lAc-7r/4)] 

- /(2A, - n/2) + - 4 - 2 j /'(2Ae - ^/2). 

This equation is difficult to solve due to the In cos term. We can use the fact that t 1 and 
expand this term about u ^2/tc — t^/2. We find 

g{u) + - ^ - ^) 9\u) ^ -| In [sin(2/tj] - ^ cot(2A,) [u ^ ' ^ 



- /(2/t, - 7r/2) + - 5 - ^) nVtc - vr/2). 

The solution is 

g{u) ^ -/(2/t,-7r/2)-|lnsin(2/g 



(A.3) 



- ^ + f ) (/'(2/tc - vr/2) + I cot(2/t,)) + 



2 TT 



(A.4) 
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where A is an integration constant we will determine later. 
At the domain wall we have 

1 a + 1 a — 1 TT /.N 

x + - = . A.6 

t t 2 ^ ^ 

The condition for the field to be constant along the domain wall becomes 

(l)jiit,x = ail/t-l/Q + l/t,-Tr/2) = 

J. / a + 1 a — 1 n\ 1 ^, f a + 1 a — 1 vr 



t tc 2 y t \ t 2 









({a -I) 




-a 









(A.7) 

2 J ' V'" \t tj 2 

— —\nt = k, 
3 

where g is given by ( ]A.4[) and k is the value of the field along the wall. This equation is 
difficult to solve in closed form, but can be solved using Mathematica. We follow a similar 
procedure as we did to obtain the solution for g. We rewrite everything so that the argument 
inside / becomes a dummy variable u. In this case, u = — — |. Then we re-express 
all t dependence in terms of u. This gives a differential equation for f{u) once we utilize 
( ]A.4p for the correct argument inside of ^ as a function of something involving u. Then one 
solves the equation for f{u). 

The general solution will have two integration constants, one each from the differential 
equations for / and g. Demanding that f{2/tc — 7^/2) is continuous at the radiation line 
fixes both constants to lowest order in 1/t: 

(p{t,x) ^ e(-x- l/t + 2/tc- vr/2)0(x-a/t- (1 -a)/tc + vr/2)0//(t,x) 

+ e(x + l/t-2/tc + vr/2)0(-x- l/t + 7r/2)0(x- l/t + 7r/2)07(t,a;),^ 

where (j)n is the solution from Mathematica and 0/ = — ^ ln(t cosx). 



B The Spherical Harmonics 

The spherical harmonics are given by 



where 



Jm 

PUx) = (-l)-(l - ^T^'^Pii^) (B.2) 
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are associated Legendre polynomials. In particular the spherical harmonics are chosen such 
that 



v* — ( —^ 

^ Im \ ^) ^ l,-rm /„ „^ 

(fnYirr,{n)Yi*^, = SwSmm'- 
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